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Abstract 

An analytical formula for the probability density function (PDF) of the velocity 
fluctuation in fully-developed turbulence is derived, non-perturbatively, by assuming 
that its underlying statistics is the one based on the generalized measures of entropy, 
the Renyi entropy or the Tsallis entropy. The parameters appearing in the PDF, 
including the index q in the generalized measures, are determined self-consistently 
with the help of the observed value [i of the intermittency exponent. The derived 
PDF explains quite well the experimentally observed density functions. 
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In this paper, we analyze the probability density function (PDF) of the veloc- 
ity fluctuations in turbulence measured in the experiments conducted by Lewis 
and Swinney [1] and by Gotoh et al. [2] in terms of a generalized statistics 
constructed on Renyi's [3] or on Tsallis' entropy [4,5]. The Renyi entropy has 
an information theoretical basis, and has the extensive characteristics as the 
usual thermodynamical entropy does. On the other hand, the Tsallis entropy 
is non-extensive, and therefore it provides us with one of the attractive trials 
for generalized statistical mechanics which deals with, for example, systems 
with long range correlations. Note that the distribution functions which give 
an extremum of each entropy have a common structure (see (5) for the present 
system) in spite of different characteristics of these entropies. The investigation 
of turbulence based on the generalized statistics was set about by the present 
authors [6] with an investigation of the p-model [7,8]. Then, the treatment 
was sublimated by a self-consistent analysis following the usual procedure of 
ensembles in statistical mechanics, where it was shown that the derived ana- 
lytical expression for the scaling exponents of the velocity structure function 
explains quite well the experimentally observed exponents [9-12]. It has been 
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revealed also that, within the same line of the statistical mechanical treat- 
ment, the PDF of the velocity fluctuations can be derived analytically [12,13]. 
An input for the ensemble theoretical method is the experimentally measured 
intermittency exponent. The parameters appearing in the ensemble theory 
are determined self-consistently with the help of the physically meaningful 
equations, i.e., the conservation of energy, the definition of the intermittency 
exponent and the scaling relation relevant to the analysis by statistics based 
on the generalized entropy [9-13]. The present ensemble theoretical method 
rests on the following aspects: 1) the fully developed turbulence is the situa- 
tion where the kinematic viscosity can be neglected compared with the effect 
of the turbulent viscosity that is a cause of intermittency; 2) in this situation, 
one can find an invariance of the Navier-Stokes equation in certain scale trans- 
formation which contains an arbitrary exponent that measures the degrees of 
singularity in velocity gradient; 3) for the distribution of the exponent, the 
statistics based on the generalized entropies is adopted, which was shown [10- 
13] to be a generalization of the statistics in the log-normal model [14-16]. The 
statistical mechanical procedure based on the Boltzmann-Gibbs entropy gives 
the results of the log-normal model, whereas the one based on the generalized 
entropies provides us with the above mentioned new results [12]. 

The basic equation describing fully developed turbulence is the Navier-Stokes 
equation du/dt + (u ■ V)u = — V (p/p) + zA7 2 m of an incompressible fluid, 
where p, p and v represent, respectively, the mass density, the pressure and 
the kinematic viscosity. Our main interest is the correlation of measured time 
series of the streamwise velocity component u, say re-component of the fluid 
velocity field u. Within Taylor's frozen flow hypothesis, the quantity of our 
interest reduces to 5u(r) = \u(x + r) —u(x)\ representing the spatial difference 
of the component u. Then, the Reynolds number Re of the system is given by 
Re = 5u £ /is = (Io/Ik)^ 3 - Here, we used Su = (e£ ) 1 ^ 3 with the Kolmogorov 
scale [17]: Ik = (^ 3 /e) 1 / 4 (= 77), and the energy input rate e to the largest 
eddies of size £ . Here, we introduced the notation 5u n = Su(£ n ) representing 
the velocity difference across a distance r ~ £ n . 

For high Reynolds number limit Re 3> 1, or for the situation where the effect 
of the term containing the kinematic viscosity in the Navier-Stokes equation 
can be effectively neglected, we see that the equation has an invariance under 
the scale transformation [8]: f — > r u — > u ', t — > t' and (p/p) — > (p/p) with 
f' = \f,u' = A a/3 M, t' = \ 1 ~ a/3 t, (p/p)' = A 2a/3 (p/p). Here, a is an arbitrary 
real quantity. Then, we have 

6u n /5u = (£ n /£ ) a/3 , e n /e = (L/^f 1 (1) 

where e n represents the rate of transfer of energy per unit mass from eddies of 
size £ n to eddies of size £ n +i- We have put eo = e. The first equation in (1) leads 
to a singularity in the velocity gradient [18] \du(x)/dx\ = lim^ n ^ \ u ( x + £ n ) — 
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u(x)\/£ n = lim^^o 8u n /£ n for a < 3, since 5u n /£^ 3 = constant. The degrees 
of the singularity are specified by the exponent a. The singularities in velocity 
gradient fill the physical space of dimension d with a fractal dimension F(a). 
Then, the probability P(a)da to find a point to be singular is specified by the 
distribution function of singularities P(a) oc 5%~ F ( a \ which leads through the 
second relation in (1) that the energy transfer rate e n is a fluctuating quantity 
in space. The conservation of energy in the energy-transfer cascade requires 
that 

(^) = 6 (2) 

where (• • •) is defined by (■ • •) = j da - ■ ■ P(a). The intermittency observed in 
the systems of fully-developed turbulence is attributed to the distribution of 
the singularities in each different size £ n = £oS n (n — 0, 1, 2, • • •) of constituting 
eddies. Here, £ n represents the diameter of the nth eddy, and is determined by 
the number of steps n of the cascade [19]. It is assumed by putting 5 n = 5~ n 
(5 = 2) that, at each step of the cascade, eddies break into 5 pieces with 1/5 
diameter. The definition of the intermittency exponent \x is given by 

(4> = (3) 

It is convenient to introduce multifractal spectrum f(a) by f(a) = F(a) — 
(d—1) which gives us P(ct) oc 8]^^ a \ 

Regarding an eddy as a cluster having the same size as the eddy [6,12], we 
can derive the scaling relation [9-12] 

1/(1 - q) = l/a_ - l/a+, (4) 

where a± are defined by f{ct±) = 0. This scaling relation is a manifestation of 
the mixing property. Note that (4) is a generalization of the scaling relation 
derived first in [20,21] to the case where the multifractal spectrum f(a) has 
negative values [22]. 

Let us find out an appropriate distribution function P( n ' (a) oc [PW (a)] n of the 
degrees of singularity a at the nth step of the cascade. Here, we are assuming 
the independence of each step. In order to discover P^(a), we will take an 
extremum of the measure of an appropriate generalized entropy with the two 
constraints, i.e., the conservation of probability: / daP(a) = constant, and the 
fixing of the g-variance: <j 2 q = ((a — ao) 2 ) q = (j daP(a) q (a — ao) 2 )/ j daP(a) q . 
By adopting the measure of the Renyi entropy [3] or of the Tsallis-Havrda- 
Charvat (THC) entropy [4,23] as an appropriate one, we have the distribution 
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function of singularities in the form [9-13] 

pW(a) da = (Z^y 1 [l - (a - a„) 2 /(4a) 2 ] n/(1 ~ ff) da (5) 

with Aa = (2X/ [(1 — q) ln2]) 1/ ' 2 and an appropriate partition function Z^ n \ 
Note that since the Renyi entropy and the THC entropy reduce to the Boltzmann- 
Gibbs entropy in the limit q — > 1, the distribution function (5) reduces to the 
Boltzmann-Gibbs distribution function in that limit. The multifractal spec- 
trum corresponding to the distribution function of singularities is given by [9- 
13] 



f(a) = 1 + (1 - q)- 1 log 2 1 - (a - a f I {Aaf (6) 



where we put /(ao) = 1, since the fractal dimension of the multifractal sets 
is one for the present experimental setup of the system. Note that the mul- 
tifractal spectrum does not depend on the number n of steps in the cascade, 
which is a manifestation of the assumption that each step is independent. 

The mass exponent r{q) is defined by 5~ d ((ejffi) = (e£$) iI 6Z r( * H ' l9 ~ 1){d ~ 1) 
where e n £ d is the total dissipation in volume £ d , and the factor 5~ d is the num- 
ber of boxes of volume £ d in the region of volume £q.Q Substituting P(a) oc 
5 l n - f{ -^ into the definition, we have 5~ T ^ = 5~ d J"™; da'P(a') S^ a '-^ oc 
5^ a? ~^ a ^ by estimating the integral in the limit 5 n < 1, where aq is obtained 
by solving 

q = df(a)/da. (7) 



Then, comparing the exponents of 5 n in the estimate of the integral, we know 
that the relation between the multifractal spectrum f(a) and the mass expo- 
nent r(q) is given by 

f(a)=qa + r(q) (8) 



with a being determined by (7). From (7) and (8), we can derive 

a = —dr(q)/dq. (9) 

Note that (7), (8) and (9) constitute a kind of Legendre transformation be- 
tween the mass exponent r(q) and the multifractal spectrum f(a). Equa- 
tion (7) is solved, analytically, to give us aq = a — 2qX/(l + JCq) with 



1 This definition is a bit different from the one in [6,9] for general d. 
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Cq = 1 + 2g 2 (l — q)Xln2. Here, we chose the solution which is finite for 
q — > 1. Then, a max = a(<? = — oo) and a min = a(g = +00) are given by 
« mffl — a = a — a min = Aa. We now know that a 2 = (Z\a) 2 (l — q)/(3 — q). 
The conditions (2) and (3) can be written, respectively, in terms of the mass 
exponent as r(l) = and /i — 1 + t(2). 

The three quantities ao, X and g are determined by the three independent 
equations (2), (3) and (4), once the value of the intermittency exponent \i is 
given. The scaling relation (4) can be solved, by substituting a± = a 
satisfying / T (a±) = 0, in the form V2X = [(a 2 + (1 - q) 2 ) 1 ' 2 - (1 - q)]/Vb, 
or a = [2bX + 2(1 - q)V2bX] 1 / 2 , with b = (1 - 2-^) /[{I - q) In 2]. 

Let us now derive the PDF for the velocity fluctuation IJ^(x n ) with —00 < 
x n < 00, where \x n \ = 5u n /5u . It may be reasonable to assume that the 
probability n^ a \x n )dx n to find the scaled velocity fluctuation x n in the range 
x n ~ x n + dx n can be divided into two parts [13]: 

n^ n \x n )dx n = n^\ Xn )dx n + nt\\x n \)dx n . (10) 



The first term II^ l \x n ), assumed to come from thermal or dissipative fluctu- 
ation caused by the kinematic viscosity, will be considered later in this paper 
after estimating the moments of the velocity fluctuation. The second term 
IIg n \\x n \) stems from the multifractal distribution of singularities in the ve- 
locity gradient given in (5), i.e., II^ n \\x n \)dx n = P( n \a)dot. This contribution 
may be related to the fluctuations caused by the turbulent viscosity. The con- 
dition J^dXnll^^Xn) = 1 determines the normalization factor here. 
For this singular part, the relation between \x n \ and a is given by (1), i.e., 
\x n \ = <5"/ 3 , and, therefore, we have an explicit expression of IJg n \\x n \) in the 
form [12,13] 

nt\\x n \)dx n = [3/ (z^\\n6 n \)]G^(\x n \)dx n (11) 



with G in \x) = aT^l - (a(x) - a ) 7(Aa) 2 ] n/(1_<?) and a(x) = 3 In xj In 8 n . 
Note that the singular part of the PDF has a dependence on \n.5u n . 

The mth moment can be calculated as [13] {(\x n \ m )) = Jf^, dx n \x n \ m IJ^ n \x n ) = 
2 7 £> + (l-2 7 S n) )a m 6$r with 2 7 W = JZodXn \x n \ m Il£\x n ) and a 3q = 
(2/[ v /Q(l + y^)]) 1 / 2 . Here, ( m = a m/3 - 2Xm 2 /[9(l + ^CQ] - [1 - 

log 2 (l + Y /, C m / 3 )]/(l — q) is the so-called scaling exponents of velocity struc- 
ture function, whose expression was derived first in [9-13]. 

In the following, we will derive the symmetric part of the PDF of the veloc- 
ity difference. The origin of the asymmetry in the PDF may be attributed to 
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Fig. 1. Experimentally measured PDF of the velocity fluctuations by Lewis and 
Swinney [1] for R\ = 262 (Re = 540 000) and fi = 0.28 are compared with the 
present theoretical results TI^ n \^ n ) given by (13) and (14). Closed circles are the 
symmetrized points obtained by taking averages of the left and the right hand sides 
of measured PDF. For the experimental data, the distances r/rj = £n/v are j from 
top to bottom: 11.6, 23.1, 46.2, 92.5 208, 399, 830 and 1440. Solid lines represent the 
curves given by the present theory with the self-consistently determined parameters 
ao = 1-162, X = 0.334 and q = 0.471. For the present theoretical curves (13) and 

(14) , the number of steps in the cascade n are, from top to bottom: 14, 13, 11, 10, 
9.0, 8.0, 7.5 and 7.0. Dotted lines are the curves given by directly adjusting the PDF 

(15) for 9 = 1, whereas dotted-dashed lines are the ones for 9 = 2 — q'. For these 
adjusted PDF's, the indexes q h s depend on the distances rjr\, and the adjusted 
values of q' are, from top to bottom: 1.168, 1.150, 1.124, 1.105, 1.084, 1.065, 1.055 
and 1.038. For better visibility, each PDF is shifted by -1 unit along the vertical 
axis. 

dissipative evolution of eddies or to experimental setup and situations. Let 
us introduce new variable £ n = x n / ((xl)) 1 ^ 2 scaled by the variance of the ve- 
locity fluctuation, and introduce the probability density function lJg n \c, n ) for 
the singular part with this new variable through the relation lJ^ n \\C, n \)d^ n = 
II^ n \\x n \)dx n . It is a symmetric function satisfying TIg n \^ n ) = flg n \— £„). 
The relation between the new variable and a is given by |£„| = 5u n / {du 2 ) 1 ^ 2 = 
£n<*" /3-&/2 with = [272 n) ^ C2 + (1 - 2^ n) )a 2 }^ 2 and the explicit form of 
flg n \\^ n \) in terms of a reduces to 

4 (n) (|e„|) = /jW^-a/S+WCa) (12) 
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with i7 s (n) = 3(1 - 2^ ll) )/(2£ n ^27iX\\n5 n \). We see from (12) that when the 
condition (2/2 — a/3 + 1 — f{a) = is satisfied, the strong n-dependence 
of i7g n ^(^ n ) may disappear. Corresponding to the two solutions a* and 
satisfying a* < ofi for this condition, we have £t = £ n <5" t//3_< > 2 / 2 and £* = 
^ 5? /3 " C2/2 satisfying < ft < f*. 




Fig. 2. Experimentally measured PDF of the velocity fluctuations by Gotoh et al. 
[2] for R\ = 381 and \x = 0.25 are compared with the present theoretical results 
n^ n \in) given by (13) and (14). Closed circles are the symmetrized points obtained 
by taking averages of the left and the right hand sides of measured PDF. For the 
experimental data, the distances r/rj = i n jr\ are, from top to bottom: 2.38, 4.76, 
9.52, 19.0, 38.1, 76.2, 152, 305, 609 and 1220. Solid lines represent the curves given 
by the present theory with the self-consistently determined parameters ao = 1.445, 
X = 0.298 and q = 0.413. For the present theoretical curves (13) and (14), the 
number of steps in the cascade n are, from top to bottom: 21, 19, 17, 14, 12, 10, 
9.0, 8.0, 7.0 and 6.0. Dotted-dashed lines are the curves given by (15) for 6 = 2 — q'. 
For this directly adjusted PDF, the indexes q n s depend on the distances r/rj, and 
the adjusted values of q' are, from top to bottom: 1.215, 1.195, 1.175, 1.145, 1.115, 
1.105, 1.095, 1.060, 1.030 and 1.005. For better visibility, each PDF is shifted by -1 
unit along the vertical axis. 

By making use of the point f* having only weak n-dependence, we divide 
into two P arts: ^ (n) (U = for |£ n | < C, and h^(Q = 

i7*< (Cn) f° r £n — l£n|- F° r smaller velocity fluctuation, we assume that the 
PDF has Gaussian form: = ft™ (£„) + fl^dU) = {Z^)-^^, 

whereas, for larger velocity fluctuation, we assume that it is given by the 
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distribution function for singularities [13]: 



(13) 



Requiring ti<J{£n) 



and n*<(£ n ) at £* to have the same value and the same 



derivative there, we can determine the variance and the normalization fac- 
tor 4s , respectively, as T n 2 = (£*) 2 /[l+3/V)] and (4s)" 1 = H^+WW 2 
Then, we see that /7<^(£ n ) becomes [13] 

nSitn) = Ut ] exp {- [1 + 3/V)l [(UCf - l] /2} . (14) 



Figs. 1 and 2 are the comparisons of the PDF fl^(^ n ), (13) and (14), ob- 
tained by the present theory with the experiments conducted by Lewis and 
Swinney [1] for R x = 262 (Re = 540 000) and fi = 0.28, and by Gotoh et al. [2] 
for R x = 381 and \x = 0.25. In both figures, solid lines represent the curves 
by the present theory, whereas closed circles are the experimental data sym- 
metrized by taking averages of the left and the right hand sides of measured 
PDF's. In Fig. 1 (Fig. 2), the distances r/rj = £ n /i] are, from top to bottom: 
11.6, 23.1, 46.2, 92.5 208, 399, 830 and 1440 (2.38, 4.76, 9.52, 19.0, 38.1, 76.2, 
152, 305, 609 and 1220) for the experimental data, on the other hand, the 
number n of steps in the cascade are, from top to bottom: 14, 13, 11, 10, 9.0, 
8.0, 7.5 and 7.0 (21, 19, 17, 14, 12, 10, 9.0, 8.0, 7.0 and 6.0) for our theoretical 
lines. We observe an excellent agreement between the measured PDF and the 
analytical formula of PDF derived by the present self-consistent theory. 

The present analysis is based on the assumptions that the distribution of the 
singularity indexes a is multifractal and that its statistics is based on the gen- 
eralized measures of entropy, i.e., the Renyi entropy or the THC entropy. Note 
that within the present analysis we cannot distinguish which of the measures 
of entropies is responsible for the system of turbulence. The latter may be one 
of the exciting future problems both for theoreticians and for experimentalists. 
It is also one of the attractive future problems to find out some connection be- 
tween the present ensemble approach and the one based on dynamical aspect 
starting with the Navier-Stokes equation, e.g., the dynamical renormalization 
group approach, the mode-mode coupling theoretical approach and so on. 

Before closing the paper, we should mention something about the usual proce- 
dures to analyze experimentally measured PDF based on generalized statistics, 
e.g., Tsallis statistics which provides us with one of the possibilities for an ab- 
stract statistical mechanics for non-Gibbsian statistics that are being observed 
in various fields (refer to [5] for more information on recent progress). Usually, 
one is apt to write down Tsallis-type PDF directly without any fundamental 
consideration when he or she is going to analyze measured PDF. Let us ex- 
plain it in the context of the present analysis of turbulence. Assuming, from 
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the beginning, that the PDF of the velocity fluctuations is of Tsallis-type, one 
may write it in the form 

n B (tn) oc {l - (1 - g')ia 2 7(5 - 3</)} V(W) (15) 



for the measured PDF symmetrized in the sense given in this paper. Here, 
q' is the index appearing in the Tsallis or the Renyi entropies which may 
have different meaning from our q used in the present paper, and 9 is an 
extra parameter determined in order to adjust (15) to the measured PDF. 
Note that the PDF (15) of the velocity fluctuations itself has a power-law 
decay for large Su n , while the present PDF does not have simple power- law 
decay but has a bit complicated decay containing In. Su n as can be seen in 
the expression (11). In Fig. 1, we put (15) with 6 = 1 by dotted lines for 
comparison which is essentially the same as done by Beck in [24] . We see that 
this PDF cannot explain the measured one. Following the analysis done in 
[25], we took 6 = 2 — q', and put (15) in Figs. 1 and 2 for this case by dotted- 
dashed lines just for comparison. As can be observed in Fig. 1, we cannot 
see big difference between our present PDF, (13) and (14), given by solid 
lines and the PDF, (15) with 9 = 2 — q', given by dotted-dashed lines, since 
the accuracy of experimental data for PDF was only of order 1CT 5 . While in 
Fig. 2 the PDF was measured up to the order of 10 -9 , and we clearly see the 
difficulty in analyzing the data by the PDF given by Beck and others in [25]. 
There may be a possibility to remedy it by adjusting the dependence of 6 on 
q' further, whose physical meaning, however, we do not know. In addition to 
this, the analysis based on (15) has another difficulty, i.e., the index q' depends 
on the distance r/rj whose dependence is listed in the captions in Figs. 1 and 
2. There might exist a possible explanation about the r/rj dependence of q' 
by relating it to the non-extensivity of Tsallis entropy [26]. In any case, the 
present investigation tells us that one should choose appropriate variables 
before applying generalized statistics. Note that the values of q' satisfy q' > 1, 
in contrast to the values of q for the present analysis where they satisfy q < 1. It 
is worth mentioning here that the present ensemble theoretical approach [9- 
13] is quite general, in that the expression of the multifractal spectrum (6) 
based on a statistical mechanical procedure may provide an example for the 
phenomenological "entropy function" S(z) with z = 1 — a in the problem of 
large deviation statistics [27]. An analysis in this direction will be reported 
elsewhere in the near future. 

The authors would like to thank Prof. C. Tsallis and Dr. A.K. Rajagopal 
for their fruitful comments with encouragement, and Prof. T. Gotoh for his 
kindness to send them the data of his numerical simulations. 
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